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Coupling an antisymmetric tensor field to the electromagnetic field in a dyonic Reissner- 
Nordström-AdS black hole background, we build a holographic model for the paramag- 
netism/ferromagnetism phase transition. In the case of zero magnetic field, the time reversal sym- 
metry is broken spontaneously and spontaneous magnetization happens in low temperatures. The 
critical exponents are in agreement with the ones from mean field theory. In the case of nonzero 
magnetic field, the model realizes the hysteresis loop of single magnetic domain and the magnetic 


susceptibility satisfies the Curie-Weiss law. 


Introduction.—The gauge/gravity duality relates a weak 
coupling gravitational theory in a (d + 1)-dimensional 
asymptotically anti de Sitter (AdS) space-time to a 
d-dimensional strong coupling conformal field theory 
(CFT) in the AdS boundary [1-4]. In recent years, many 
fascinating condensed matter phenomena have been 
found using the duality, including holographic superflu- 
ids (superconductors) [5, 6], (non-)Fermi liquids [7-9], 
Josephson junctions [10-12], superconducting quantum 
interference device [13], magnetic properties in supercon- 
ductor [14-17], and quantum phase transitions [18]. 

Magnetism is ubiquitous in many strongly correlated 
electronic systems such as high temperature supercon- 
ductors and heavy fermion metals, and plays a central 
role in quantum phase transitions in these materials. The 
gauge/gravity duality provides new approaches and per- 
spective to understanding these challenging systems. Yet 
in holographic contexts, due to various technical chal- 
lenges, models of magnetism are scarce and not exten- 
sively explored (see e.g. [18]). In this paper we propose 
a simple holographic model for magnetism, and present 
an example of the paramagnetism/ferromagnetism phase 
transition. We believe the model proposed here should 
serve as a good starting to explore more complicated 
magnetic phenomena and quantum phase transitions. 

The modelThe paramagnetism/ferromagnetism 
phase transition breaks the time reversal symmetry 
spontaneously in low temperatures (if spatial dimension 
is more than 2, it also breaks spatial rotation symmetry). 
So we need a real field to describe this symmetry break- 
ing spontaneously. Note that in the mean field theory 
of paramagnetism/ferromagnetism phase transition, the 
magnetization as the order parameter is a pseudovector 
and is coupled to the external magnetic field; and in the 
weak magnetic field, the magnetization is proportional 
to the magnetic field in magnetic materials. Therefore 
in the spirit of AdS/CFT correspondence it is natural 
to introduce a real antisymmetric tensor in bulk, its 
space-space component describes the dual field of 
magnetization. Thus we consider a simple model with 


action as 
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Here 2k? = 167G and G is the Newtonian gravitational 
constant, À and J are two constants, m is the mass of the 
real tensor field M,,,, A, is the gauge potential of U(1) 
gauge field and j” is the U(1) current coming from some 
field carrying U(1) charge. We here introduce j” into the 
action for the convenience of discussion. In fact we will 
set 74 = 0 in this paper. The antisymmetric tensor field 
M „v is the effective polarization tensor of the U(1) gauge 
field strength Fv, whose physical meaning will be ex- 
plained shortly. V(M,,,) describes the self-interaction of 
the polarization tensor, which should be expanded as the 
even power of Mav. The lowest order has two linear in- 
dependent forms. We take V (Mpv) = M,” M,” M,° Ms" 
for simplicity. The probe limit corresponds to À —> 0. 
The equations of motion for matter fields are 
Ve Py = x VEM, ph 
w= Zi uv + Jo; (3) 
V? Muy — ml Muy — JM,’ Ms" Mz, — Fup = 0. 


From the equations for gauge field, one can see that 
the polarization field can contribute U(1) current and 
charge density. If introduces an auxiliary field such 
as Hyy = Fav — A*Myv/4, the U(1) current is only asso- 
ciated with the field H,,. This is very similar to what 
happens for electrodynamics in medium, where the aux- 
iliary fields D and H are introduced corresponding to E 
and B, respectively. 

From equations (3) one can see that j” is still a con- 
served current and associated conserved charge is just the 


U(1) charge. In the planar symmetry case, its density can 
be expressed as 


Q= I *(Fuv — X’ Mpu /4) = lim r°(Fir — X°’ Mer /4). 
r= oo Too 

(4) 

In this paper, we will work in the probe limit and the 

background is a planar dyonic Reissner-Nordstr6m-AdS 

black hole. The background geometry and gauge poten- 
tial take the following forms [19] 
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where yp is the chemical potential, B is the magnetic field, 
and the AdS radius L has been set to be unity and the 
black hole horizon r} = 1. Then the temperature of the 
black hole is 


T = (3— p? — B?) /4r. (6) 


This bulk magnetic field can be regarded as the external 
magnetic field in the dual boundary theory. In this pa- 
per we will work in grand canonical ensemble, where the 
chemical potential will be fixed as a nonzero value, while 
the magnetic field can be taken zero or not. We consider 
a self-consistent ansatz for the polarization field as 


Mpv = —p(r)dt A dr + p(r)dx ^ dy. (7) 


Then the non-trivial equations for the polarization field 
in the background (5) read 
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where a prime denotes the derivative with respect to r. 
It is interesting to see that these two equations decouple 
from each other in this case. At the horizon, we impose 
the regularity condition, which implies the following re- 
lation 
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Thus once given the initial values of p and p at the hori- 
zon, one can integrate equations (8) to get the solution. 

The asymptotic behavior and free energy.— In order to 
keep the background to be asymptotic AdS, one has to 
require suitable asymptotic behaviors for p and p near 


the AdS boundary. There are two different situations in 
terms of the value of m?. When m? = —4, a logarithmic 
term appears in the asymptotic expansion of the solution 
as r — oo. We will not consider this case in this paper. 
When m? 4 —4, we have the following asymptotic be- 
haviors 
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The Breitenlohner-Freedman (BF) bound requires m? > 
—17/4 [20, 21]. Note that there is an other restriction 
on m? coming from the U(1) charge density in (4). A 
well defined probe limit demands that the charge density 
contributed by the polarization should be finite and the 
leading order is \? order. In order to check this, we have 
to treat à? as a samll quantity and solve the full equa- 
tions of motion of (1) order by order. However, since the 
charge density (4) is calculated at r — co, the asymp- 
totic AdS boundary guarantees that we need not to solve 
metric if we only care for the order \?. A direct calcula- 
tion gives the contribution coming from the polarization 
in (4) is 
pà? 


X _ 
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We see that if p+ #0 or ô < 1, the quantity is divergent, 
which leads the system to be inconsistent. In order to 
avoid this, we further impose 


p+ =0, m? > —4. (12) 


On the other hand, in order to the condensate happens 
spontaneously when the temperature is lowered, the m? 
of the real tensor field should violate the BF bound of 
the field in AdS_ x R?. This leads to the constraint of 
m? as 


3 
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Note that here AdS x R? is the near horizon geometry 
for an extremal black hole (5) with a vanishing temper- 
ature. Therefore the real tensor field M does not violate 
its BF bound in AdS4, but does violate in AdS x R?. 
As a result, the spontaneous condensate is an effect of 
temperature, it is not sensitive to the details of the self- 
interaction terms for the tensor field M. 

The free energy of the dual boundary theory can be 
obtained by calculating the on shell action of bulk in 
the Euclidean sector. However, for a generality and a 
comparison with the Ginzburg-Landau theory, we here 
present the off-shell free energy of the dual field theory. 
In other words, we calculate the bulk action when the 


equations of motion are no longer required to hold. In 
the probe limit, a direct calculation gives the free energy 
density of the polarization tensor as 
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where we have set 2x? = 1 and the operator L has been 
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here the contribution coming from the component p to 
the free energy did not included because it is not relevant 
to the following discussions. In the on shell case with the 
equation of motion, Lp = Jp?/r®—B/r?, Gonsheu diverges 
if ps #0. This divergency coming from the boundary 
terms can be canceled by adding some suitable surface 
counter term, while the one from the volume integration 
part cannot be renormalized in the case of B £ 0. So in 
order to make the on shell action finite, we need 
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p+ =0. (15) 


In this case, the boundary terms appearing in the action 
vanish. The condition (15) can also be understood as fol- 
lows. In the AdS/CFT correspondence, when B = 0, p+ 
and p_ correspond to the source and vacuum expectation 
value of dual operator in the boundary field theory (up 
to a normalization factor), respectively. Therefore one 
should take p4 = 0 since one wants the condensation to 
happen spontaneously. 

In the on shell case, the free energy density for the 
polarization field can be written as 


Nem ee Jp’ 
onshell = -F ae — BN, 1 
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where 
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According to this free energy density, the term propor- 
tional to J can be interpreted as the self-interacting en- 
ergy of polarization, and N can be explained as the mag- 
netic moment of the dual boundary field theory. 

The spontaneous magnetization in low temperature 
corresponds to the spontaneous condensate of p in the 
bulk in the absence of external magnetic field. One can 
see from (16) that once the spontaneous magnetization 
phase appears it is always more favored than the triv- 
ial solution without the condensation if J is negative. 
On the other hand, from the first equation of (8), the 
term proportional to J likes the —\¢* term in the Higgs 
mechanism, one therefore requires J < 0 in order the 
symmetry to be broken spontaneously. 


To calculate the off shell free energy, let us consider 
the following Sturm-Liouville eigenvalue problem 


Lpn(r) = Anw(r)pn(r); (18) 


with boundary conditions: p,,(7+) is finite at horizon and 
Pn(co) = 0 at r = co. Here An’s are eigenvalues and the 
weight function w(r) can be an arbitrary positive contin- 
uous function in the region r € [r;,0o) so that the eigen- 
value àn will not influence the asymptotical form of the 
solution pn at the infinity and pn is square integrable in 
the region r € [r4}, 00). For example, we may take w(r) = 
1/r* with k > 2 for simplicity in our case. Thus any func- 
tion, h(r), which is finite and square integrable in the 
r € [r4,00), can be expanded as A(r) = Xp] cnpn(r), 
where ¢, = (Pash) = I w(r)pn(r)h(r)dr. Applying 
this method to our case, we can show that near the crit- 
ical temperature Te, the free energy density can be ex- 
pressed as 


Gortshen /A? = a1 (T/T. — 1)N? + agJN* — BN + O(N), 
(19) 
where the coefficients a, and ag, and the critical temper- 
ature Te can be numerically determined by solving the 
equation Lp,(r) = 0. The details will be presented in 
[22]. One can see that the free energy density is quite 
similar to the one of the Ginzburg-Landau theory for the 
paramagnetism/ferromagnetism phase transition. 1 

The spontaneous magnetization and the response to the 
external magnetic field In order to study the behavior 
of spontaneous magnetization and the response to the 
external magnetic field B, we have to solve the equa- 
tions (8) numerically. Since our focus is on the behavior 
of magnetic moment near the critical temperature, we 
pay attention on the equation of p. As an example, we 
take J = —1, m? = —3 in equations (8). The other values 
for J and m? satisfying the restriction (13) and J < 0 
give the qualitatively same results. 

Firstly, we consider the spontaneous magnetization of 
the system. For this, we set B = 0. Solving (8) with 
requirement (15), we find there exists a critical temper- 
ature Te. When T < Te, the nontrivial solution p 4 0 
appears. The magnetic moment N as a function of tem- 
perature is plotted in Figure 1. In this figure we also 


1 According to the AdS/CFT dictionary, one should view p— as 
the magnetic moment in the dual field theory. From the form 
(19) of free energy, we think it is more natural to view N rather 
than p— as the magnetic moment in this paper since it is N which 
is coupled to the external field B rather than p_. In fact we can 
also express N as a function of p_. In the symmetric phase with 
B = 0, both N and p— are zero, and in the condensed phase close 
to the critical temperature, p— and N are indeed proportional to 
each other with the same critical behavior. Far from the critical 
temperature, the nonlinear relation between N and p— is due to 
the self-interaction term V(M“”) of the tenor field in (3). 
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FIG. 1. The magnetic moment as a function of temperature. 
Here the critical temperature is T./ œ 0.00915. 


plot the behavior of p_. Indeed we can see that near the 
critical point, N and p_ have the same behavior. 

By fitting this curve near the critical temperature, we 
find that there is a square root behavior for the magnetic 
moment versus temperature, which is a typical behavior 
for a second order phase transition. Specifically, for J = 
—1, m? = —3 we have 


N?/M yu? ~ 6.5054 x 1073(1 — T/T). (20) 


This gives the critical exponent 1/2, the same as the one 
from mean field theory. Note that if we take p_ as the 
order parameter, it has the same critical behavior near 
the critical temperature. The second order phase transi- 
tion can also be checked by computing the free energy 
of the system. It is easy to see that the free energy 
is continuous at the critical temperature. Thus we see 
that when T < Te the magnetic moment occurs spon- 
taneously, which leads to the breaking of time reversal 
symmetry spontaneously. It is natural to expect that 
this system enters into a ferromagnetic phase as T < Te. 

Next let us turn on the external magnetic field B to 
examine the response of magnetic moment N. This can 
be described by magnetic susceptibility density x, defined 
as 


N 
x= lm =] . (21) 


In this model, the magnetic susceptibility comes from 
two parts, one is from the background black hole and 
the other from the polarization field. The pure dyonic 
Reissner-Nordström-AdS black hole corresponds to a dia- 
magnetic material [23]. Although in the probe limit, 
the contribution from the polarization field is suppressed 
by A?, near the critical temperature, we find that the part 
from the polarization field dominates for any nonzero À? 
as this part is divergent at the critical temperature. In 
the right plot of Figure 2, we show the magnetic suscep- 
tibility as a function of temperature and find it satisfies 
the Curie-Weiss law of ferromagnetism in the region of 
T> Ts 


xd? /u œ 0.5463(T/T, — 1). (22) 
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FIG. 2. Left: The magnetic moment with respect to external 
magnetic field B in the different temperature. Right: The 
magnetic susceptibility as a function of temperature. 


Now we can conclude that the dual system is in a param- 
agnetic phase in high temperatures and a ferromagnetic 
phase in low temperatures. The model (1) can describe 
a paramagnetism/ferromagnetism phase transition. 


In the left plot of Figure 2, we show the magnetic mo- 
ment with respect to external magnetic field B in the 
different temperatures. In the case of T < Te, we can 
see that the magnetic moment is not single valued. The 
parts DE and BA are stable, which can be realized in 
the external field. The part CF is unstable which cannot 
exist in the realistic system. The parts EF and CB are 
metastable states, which may exist in some intermedi- 
ate processes and can be observed in experiment. When 
the external field continuously changes between — Bmax 
and Bmax periodically, the metastable states of magnetic 
moment can appear. Thus we see a hysteresis loop A-B- 
C-D-E-F-A. This corresponds to the hysteresis loop in 
the single magnetic domain. 


Summary and Discussion._We have built a simple 
(3 + 1)-dimensional bulk gravity model which can re- 
produce some salient properties of a (2 + 1)-dimensional 
paramagnetism/ferromagnetism phase transition. In the 
case of zero external magnetic field, the model can break 
the time reversal symmetry spontaneously in low temper- 
atures and realize the paramagnetism/ferromagnetism 
phase transition. The critical exponent agrees with the 
one from mean field theory. In the case of nonzero mag- 
netic field, the model realizes the hysteresis loop of single 
magnetic domain and the magnetic susceptibility satisfies 
the Curie-Weiss law. Furthermore, let us mention here 
that the relation (20) for the magnetic momentum and 
(22) for the magnetic susceptibility can also be obtained 
through the off shell free energy (19) [22]. 

Here some remarks are in order. In the ansatz (7) we 
have not considered the components Mj; and Mty. In 
fact, these two components can be viewed as dual field 
of the order parameter in the paraelectric/ferroelectric 
phase transition in dielectric materials. The spontaneous 
electric polarization can be realized in this model, which 
will be studied in a forthcoming paper [22]. As a re- 
sult, our model provides a unified holographic descrip- 
tion for the paramagnetism/ferromagnetism and para- 
electric/ferroelectric phase transitions. In addition, we 


can extend this model by including two different polar- 
ization tensors in order to form two different magnetic- 
ordered structures. By adding some suitable interaction 
between them, one can expect that the antiferromag- 
netism and ferrimagnetism phases will appear [24]. In 
particular, we would like to stress that this paper is just 
a first step to stress magnetic orders in some strongly 
correlated systems. By combining this model with some 
holographic superconductor models, one can explore the 
competition and coexistence between magnetic-ordered 
phase and superconducting phase in high temperature 
superconducting materials. 
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